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Non-minimal matter couplings have recently been considered in the context of massive gravity 
and multi-gravity. These couplings are free of the Boulware-Deser ghost in the decoupling limit 
and can thus be considered within an Effective Field Theory setup. Beyond the decoupling limit 
the ghost was shown to reemerge in the metric formulation of the theory. Recently it was argued 
that this pathology is absent when formulated in terms of unconstrained vielbeins. We investigate 
this possibility and show that the Boulware-Deser ghost is always present beyond the decoupling 
limit in any dimension larger than two. We also show that the metric and vielbein formulations 
have an identical ghost-free decoupling limit. Finally we extend these arguments to more generic 
multi-gravity theories and argue that for any dimension larger than two a ghost is also present in 
the vielbein formulation whenever the symmetric vielbein condition is spoiled and the equivalence 
with the metric formulation is lost. 


Introduction. - With the successful development of 
ghost-free massive gravity [1, 2], bi-gravity [3] and 
multi-gravity [4] theories (see Refs. [5, 6] for reviews), 
the question has naturally arisen of whether there are 
more general consistent ghost-free couplings to matter. 
Unlike GR, non-minimal couplings to both metrics are 
not forbidden by any symmetry, and will inevitably arise 
in any low energy effective field theory of massive gravity 
or multi-gravity. A number of possible phenomenological 
couplings to matter have been explored in the literature 
[7-20] and their theoretical consistency was explored in 
[21-28], 

As is well known, all multi-gravity theories have a 
low strong coupling scale A which lies parametrically 
between the graviton mass and the Planck mass. For 
massive gravity theories with a fixed reference metric, 
in D dimensions the scale is A = ( m 2 M^ 2 )/ 2 ) 2 /(-°+ 2 )_ 
For multi-gravity it is a generalization which is typically 
dominated by the lowest non-zero graviton mass and the 
lowest Planck scale [29, 30]. If these theories are viewed 
in terms of an effective field theory (EFT), then A may 
be taken to be the cutoff of the EFT. In which case, we 
only require for consistency of the theory, the absence 
of ghosts below the cutoff scale A. In particular any 
non-minimal matter coupling which remains ghost-free 
in the decoupling limit (DL) Mp\ oo keeping A fixed 
is considered acceptable from an EFT point of view. In 
a recent work precisely such a coupling was proposed in 
the metric formulation of massive gravity and bi-metric 
theories, and amounts to an effective ‘composite’ metric 
g e g in [22], 

dfv = 9 V + 2a 9n» (vV 1 /) ^ + a 2 f^ , (1) 

where and are two metrics (in bi-gravity they 
are the two dynamical metrics, while in massive grav¬ 
ity 9p.v can represent the dynamical metric and the 
Minkowski reference metric). 

In terms of the effective metric (1), the coupling to 


matter then takes the standard form, 

•^matter — £g {9,dlpg,1pg)+Cf (/, dtpf^f) 

+ Coff {g e s,dx,X ) , (2) 

where ip g , iff and x symbolize different matter sectors. 
The generalization of the coupling to the effective met¬ 
ric (1) was also considered for multi gravity in [9]. In 
the vielbein formalism 1 the effective vielbein takes the 
remarkable simple form 

N 

e e ff = ^a/e (/) , (3) 

i=i 

where e^ is the vielbein for the metric g$, and the ct/ 
are arbitrary dimensionless coefficients. 

Vainshtein. - As we have emphasized, from an EFT 
point of view the absence of ghosts below A is a suffi¬ 
cient criterion. Nevertheless, it has been traditional to 
ask for a tighter restriction in massive gravity theories. 
If the additional interactions were ghost-free nonlinearly, 
then the Vainshtein mechanism allows us to make sense 
of the theory above the scale A, by redressing the cut¬ 
off scale around certain backgrounds. Thus in order to 
address the viability of the Vainshtein mechanism, it is 
natural to ask whether there exist any non minimal cou¬ 
plings to matter which are ghost-free nonlinearly. In the 
case of the composite metric proposal of [22] this was 
definitively shown to be ghostly [24, 25]. More recently, 
it has been argued that the unconstrained vielbein ver¬ 
sion of this interaction was ghost free [27]. A priori this 
is not a contradiction since the vielbein and metric for¬ 
mulations are not equivalent. This is because integrating 
out Lorentz transformations that encode the extra vari¬ 
ables in the vielbein no longer imposes the symmetric 


1 See Refs. [4, 31, 32] for the formulation of massive gravity in the 

vielbein. 
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vielbein (or Deser-van Nieuwenhuizen) condition [33, 34], 
but a non-trivial condition which depends on the matter 
content. As such the metric formulation of the uncon¬ 
strained vielbein formulation would contain highly non¬ 
trivial non-minimal matter couplings and would differ 
from a effective coupling to (1). 

Boulware Deser ghost. — Before moving to the core 
of the argument, we first clarify what we mean by the 
Boulware-Deser ghost as first discovered in [35]. The 
Boulware-Deser (BD) ghost is usually associated with 
the loss of the primary constraint generated by the lapse 
N. However one can always push the problem to the 
level of the secondary constraints by introducing auxil¬ 
iary variables. For instance one may consider the Hamil¬ 
tonian H = NRq+N 2 Q , where Rq and Q are functions of 
the dynamical phase space variables. In this formulation 
it is clear that N does not generate a primary constraint. 
For this example, we can restore that primary constraint 
at the price of spoiling the secondary constraints by in¬ 
troducing an auxiliary variable a and conjugate momen¬ 
tum p„, with H = N(Rq + crQ) — \cr 2 Q + A p a , where 
A is a Lagrange multiplier that ensures that p a = 0. In 
that formulation, the primary constraint associated with 
the lapse N survives C\j = Rq + crQ ss 0. In addi¬ 
tion, there we have another primary constraint namely 
C‘ 2 .[ = p a ~ 0. The problem manifests itself at the level 
of the secondary constraints, where the one associated 
with C 2 j now involves the lapse: 

C 2 ,n = {C 2tI ,H} = -— =q(N--o\ «0, (4) 

and this equation can be solved for the lapse N rather 
than one of the dynamical variables and so is no longer 
a true constraint. Similarly the secondary constraint 
associated with C\j can now be solved for the Lagrange 
multiplier A and we thus loose two secondary constraints. 

In the formulation of [27] the primary constraint asso¬ 
ciated with N exists in all dimensions, but the existence 
of secondary constraints needed to ensure the correct 
number of degrees of freedom in D > 2 dimensions 
are not guaranteed 2 . In that case the analogue of the 
auxiliary variable a are the Lorentz rotations (to be 
defined below) and the number of secondary constraints 
associated with them increases with the number of 
dimensions. In this language we may thus wonder if 
the number of ghost would increase with the number 
of dimensions. The answer is no. The reason for that 
is while all the L(Z? — 2)(D — 1) secondary constraints 
associated with the rotations involve the lapse, it is 


2 We thank K. Hinterbichler and R. Rosen for discussions on this 
point. 


always possible to reformulate them in a way where 
only one of these equations depend on the lapse and 
the others are free of the lapse and therefore genuinely 
generate ^(D — 2) (D — 1) — 1 secondary constraints. The 
second secondary constraint which is lost is the same 
one as that associated with C\j which can be solved 
for one the Lagrange multipliers A and so a total of 
two secondary constraints are lost. So even though the 
number of auxiliary variables increases with the number 
of dimensions, the number of ghosts does not. In this 
sense the ghost found here can be attributed to the BD 
ghost as opposed to the ghosts found in [36] and in [37] 
which do increase with the number of dimensions. 

In this manuscript we shall diagnoze the problem di¬ 
rectly at the level of the primary constraint by integrat¬ 
ing out the auxiliary variables (namely the rotations). In 
this language it is more manifest that we are still dealing 
with the BD ghost, i.e. a failure of the Hamiltonian to 
be linear in N. 

Counting. - To understand the loss of the BD con¬ 
straint, let us first consider the massive gravity limit. 
Since massive gravity arises as a consistent decoupling 
limit of multi-gravity theories, the existence of a ghost 
in the massive gravity limit implies a ghost in any multi 
gravity extension by the usual reasoning. 

In what follows we use the notation that greek in¬ 
dices represent space-time indices g,u = 0, • • • ,D— 1, 
i,j = 1, • ■ • ,D — 1 run over space dimensions, A, B,- ■ • = 
0, • • • , D — 1 run over all the Lorentz directions while 
a, b, ■ ■ ■ = 1, • • • , D — 1 run over only the space Lorentz 
directions. 

We begin with the vielbein formulation of massive 
gravity with the usual EH kinetic term and with mat¬ 
ter coupling to the composite vielbein e + af. Working 
in unitary gauge f^ A = 8 A , we may decompose the dy¬ 
namical vielbein e a into a Lorentz boost/rotation of an 
upper-triangular vielbein which contains the same num¬ 
ber of variables as the metric [4] 

e/ = A A B E I1 B (5) 

where E A is upper-triangular along the time-direction 
and symmetric in the space-space directions. Specifically 
it takes the form 

E° = Ndt (6) 

E a = N i E i a dt + E i a dx i , (7) 

where £)“ = E *. The spatial 3-metric is gij = 
E^EjSab- In general we may factorize the Lorentz 
transformations into 

A a b = A 2 ^ (boost) A C B (rotations), (8) 

where the rotations act only in the space-space direc¬ 
tions. 
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For instance in D = 3 dimensions, we express the nine 
components of the vielbein in terms of a lapse TV, two 
shift TV*, two boosts v a , one rotation r and the three 
components of a symmetric spatial vielbein E i a = 

A _ ( TV 7 + N i e i a v a Nv a + tfe* + ^v b v a ^j \ 

^ ~ { e“v a e? v b v a ) j ’ 

where 7 is given by 7 = \/l + v 2 and e° is the rotated 
version of E i a , 


= (A r ) a b E t 


7r r ) ( E L 3 E* \ 

-r lr )\ E ! 2 E 2 2 ) 


( 9 ) 


with 7 r = v 1 — r 2 . In arbitrary dimensions, this decom¬ 
position may be summarized as 

D 2 (vielbein) = (D — 1) (boosts) (10) 

+ —(D — 2)(D — 1) (rotations) 

+ D (lapse and shift) 

+ I-D(D — 1) (spatial metric), 

where the spatial metric is given uniquely in terms of 
the spatial symmetric vielbein £7“. 


Now the crucial property is that as long as we choose 
the EH kinetic term, the (D — 1) boosts and the — 
2)(D — 1) rotations drop out of the kinetic term, mean¬ 
ing that in unitary gauge they are non-dynamical vari¬ 
ables which do not enter the phase space symplectic form. 
In addition, the lapse and shift TV, TV * are also non- 
dynamical, leading to a total of ^D(D + 1 ) auxiliary 
variables. The \D{D — 1) E i “ are dynamical and after 
integrating out the spin-connection all receive a kinetic 
term. If none of the auxiliary variable played the role of 
a Lagrange multiplier, the total dimension of the phase 
space would be 

Phase space dimension = 2x ^—D(D — 1)^ . (11) 

On the other hand the number of degrees of freedom of 
a massive graviton in D dimensions, is that of a massless 
one in D + 1 dimensions, namely \D(D — 1) — 1. 
As a result we would obtain in any dimension one 
additional degree of freedom which is the BD ghost. We 
thus require for the absence of the BD ghost that one 
combination of the auxiliary variables acts as a Lagrange 
multiplier. 


In [27] the existence of this additional constraint was 
argued as follows. Since the Hamiltonian is clearly lin¬ 
ear in the shift TV 1 , we may use the shift constraint to 
solve for the D — 1 Lorentz boosts in a manner which 


preserves linearity in TV. This then ensures that TV im¬ 
poses a constraint, and so naively the BD ghost must be 
absent. This argument is correct in D = 2 since their are 
no rotations. However in all dimensions D > 2, it is nec¬ 
essary to further check that it is possible to integrate out 
the rotations in a manner which preserves the linearity 
in TV. Equivalently we may use the fact that since the 
rotations do not have a kinetic term, we can introduce 
a momentum conjugate for them and an additional pri¬ 
mary constraint that they have vanishing momenta. For 
the counting to work out, these must lead to secondary 
constraints. The secondary constraints associated with 
the vanishing of the momentum conjugate for the rota¬ 
tions are none other than the space-space part of the 
symmetric vielbein condition and the requirement that 
these are truly constraints is equivalent to asking that 
the rotations can be integrated out in a manner which 
preserves linearity in TV. The arguments of [27] assumed 
the secondary constraints which fix the rotations would 
arise and this is what we find not to be the case. 

Symmetric Vielbein Condition. In the case of mas¬ 
sive gravity with minimal coupling to matter the sec¬ 
ondary constraint exists and this can be understood as 
follows: As is well known, in the case of minimal cou¬ 
pling to matter, varying the action with respect to the 
boosts and rotations reproduces the symmetric vielbein 
condition 3 


e^f v B r]AB = e v A S B r]AB ■ ( 12 ) 

If we consider only the space-space part of this equa¬ 
tion, e i A fj B r)AB = e A fi B VAB then those constitute 
|(D — 2)(D — 1) equations. It is precisely these space- 
space parts that arise as the secondary constraints 
associated with the vanishing of the momentum conju¬ 
gate to the rotations. 

Crucially since e A is independent of the lapse and 
shift, then these \(D — 2){D — 1) equations may be 
solved algebraically to determine the A(AT> — 2 )(ZD — 1 ) 
Lorentz rotations in terms of E i “ and the boosts v a , 
and these solutions are independent of TV and TV*. In 
other words these are truly constraints on the system 
which can be used to fix the rotations. Substituting 
back into the action, the Hamiltonian remains linear in 
TV 1 and the associated TV* constraints may be used to 
solve for v a . This leaves behind the TV which imposes 
the final primary constraint which removes the BD ghost. 


3 There generically exist other branches of solutions of these 
equations in which the symmetric vielbein condition does not 
hold (see for example [38, 39]), however these are disconnected 
branches in which ghosts are present. 
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In the case of the non-minimal matter coupling, the 
symmetric vielbein condition is lost and gets sourced in¬ 
stead by the stress-energy tensor of the non-minimally 
coupled matter fields. Since this stress-energy depends 
on the shift and lapse, it follows that these ^(D—2)(D—1) 
equations can be solved to determine the ^(D — 2)(D—1) 
Lorentz rotations in terms of not only a and the boosts 
v a , but also N and N l . This dependence on the La¬ 
grange multipliers means that these are not true con¬ 
straints any more (or at least not all of these equations 
are constraints). Plugging these expressions back in the 
Hamiltonian, the latter will no longer be linear in the 
shift and the lapse, hence spoiling the usual constraint 
that removes the BD ghost. This argument works in all 
dimensions D > 2 and for arbitrary matter fields. For 
concreteness, we focus in D = 3 dimensions, with the 
understanding that adding dimensions to the theory will 
not alleviate the problem (quite the opposite). 

Integrating out the auxiliary variables in massive 
gravity. In order to better compare with the case of 
the non-minimal coupling, we start by looking at the 
mass term explicitly. To simplify the discussion, we con¬ 
sider the following mass term 

2 

Hmass = Mpi-^-eABC e A A f B A f c . (13) 

Expressed in terms of the dynamical and auxiliary vari¬ 
ables, the Hamiltonian for this mass term takes the form 


H mG R = N°R + N l Ri 


+Mp\m 2 


Nl + ef [% + 


1 + 7 


+ N'efv* 


(.14) 


We can now integrate out the auxiliary variables as fol¬ 
lows: We first solve the equation of motion for the shift 
in terms of the boost v l which will then be expressed in 
terms of the dynamical variables E { “ and the rotations, 

Mpim 2 e i a v a = -Ri. (15) 


Even though Vj depends on the rotations, the rotation 
matrix A r defined in (9) drops out of the inner product 
v 2 . As a result the Lorentz factor 7 is independent of 
the rotation which only enter in e. ^<5® + 77 ( 7 ) ■ Since 

the lapse does not enter that expression, we can easily 
integrate out the rotation without involving the lapse. 
As a result the Hamiltonian remains linear in the lapse 
after integration of the auxiliary variables. 

Non-minimal coupling. This result only holds 
for a minimal matter coupling which does not include 
the Lorentz Stiickelberg fields. When matter is chosen 
to couple to the composite vielbein, e e g = e + af , 
this coupling breaks local Lorentz invariance, then the 
matter Lagrangian contains explicit dependence on the 
Lorentz boosts and rotations. Varying the action with 


respect to these Lorentz boosts/rotations gives now a 
modification of the symmetric vielbein condition which 
includes a contribution from the matter Lagrangian. 
However generically this condition depends explicitly on 
N and N l . We may continue to solve the space-space 
component of this equation for the Lorentz rotations, 
but the solution will now depend on N. Substituting 
back in the action, this will generate non linearities in 
N not present for minimal couplings. These destroy the 
existence of the BD constraint as we shall see explicitly 
below. Rather than re-deriving the modification of 
the symmetric vielbein condition, we follow the same 
procedure as previously and first solve the shift equation 
for the boosts and then integrate out the rotation. 


For definiteness we consider a scalar field 7 . Its po¬ 
tential goes as yj—g e g V (%) and can always be absorbed 
in a redefinition of the mass terms and is thus ignored 
here. Moreover to go straight to the point, we focus on 
the ultra local limit. This is justified since the gradient 
terms for matter cannot possibly compensate the contri¬ 
butions from the kinetic term. In other words, if a con¬ 
straint is lost in the ultra-local limit then it is lost in the 
whole theory. Then in terms of the momentum conjugate 
p x associated with the scalar field x, the contribution to 
the Hamiltonian from the non-minimal matter coupling 
is 


Hit 


1 

2 


N e S 2 

, P x , 


(16) 


where N e g is the effective lapse of the composite vielbein 
< 7 ®® is the effective spatial metric. These effective quan¬ 
tities can be expressed in terms of the shift, lapse, boost, 
rotation and spatial symmetric vielbein of the dynamical 
vielbein (see for instance Eq. (7) of Ref. [27]). As men¬ 
tioned in Ref. [27], the effective lapse N e s is linear in the 
shift and lapse N l and N and gf® is independent of the 
shift and lapse. Notice however that both A r e ff and gf® 
depend non-trivially on the rotations. 

The exact expression of the effective quantities in terms 
of the dynamical metric is complicated but for the pur¬ 
pose of this discussion, it is sufficient to perform an ex¬ 
pansion in a (recalling that a = 0 corresponds to the 
minimal matter coupling). Once again for concreteness 
we also work in D = 3 dimensions, although the results 
are generalizable to any D > 2 dimensions. 

In terms of the lapse, shift, boost and spatial vielbein, 
the non-minimal matter Hamiltonian (16) takes the form 
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to second order in the parameter a 

_2 „ 


H, 


1 N 


- -p\ + ^ (N' Vi + 7 ) - -- G(e, E) 

2 Vff V9 L 2 ^ 


matter — 


TV 


+ 


ct 2 p 2 


(1 - 2 7 2 )e/ + (1 + 2 7 ) ei ° 


1 + 7 


+ N‘v°eM ( — 2'y5ijS Kt + 2Si 


M 


V V 


1+7 


NQ(e, v i ) 


+ O (a 3 ) , 


(17) 


where g denotes the determinant of the spatial metric, 
g = det gtj and is invariant under rotations. In D = 3 
dimensions, g = (e y e a be i a e J b ) 2 = (e u tahE^E^) 2 . G and 
Q are functions of the boost and the rotation given by 


G = e^e ab 5“ + 


1 + 7 


(18) 


and 


Q — &ij &kt 


- (1 + 7)(8 - 37 - 6 7 2 + 2j 3 )S ij S u 


+ (12 + 7 7 - \ 1 2 )^v k v l + (1 + 7 ) (-3 + 7 )<fW 


- (1 + 7 )5 u S kj - (2 + 7 )6 ie v k v j - vVt> V 


• (19) 


Integrating out the auxiliary variables. - We may 

now integrate out the auxiliary variables as we did in the 
case of massive gravity. The shift equation leads to the 
equation for the boosts: 

p 2 

Mpim^e^Va + a—pVi ( 20 ) 

V9 

W^iPeu (s?6j - 2 7 S i0 5 u + 26 zj ^-) = -R ,, 
2 9 \ 1 + 7/ 


which can be easily solved perturbatively in a to obtain 
an expression for the boosts in terms of the symmetric 
vielbein Ejj and the rotation. Next we derive the equa¬ 
tion of motion for the rotation. To simplify this deriva¬ 
tion, we perform a perturbation in Ri. To quadratic order 
in Ri , the rotation r defined in (9) satisfies 


EUi,E k R k 

8M 2 l m±[E]Jg 


1 + 


ap~ 


2 2 
ocV x 


2 MpimV 

p 2 

8Mpim 2 V 2 


Vd 


AMpim^g 


(4iV-[^D (21) 


fl[E 2 ]-^[E] 2 +6N[E] 


l[E 2 ] + l[E} 2 -8[E}N + AN- 



+ O (/Jf/Mpjm 6 , a 3 ) 


where we used the notation [ E] = £+ and [E 2 ] = ElEf. 
We now see explicitly that the rotation depends on the 
lapse N and so integrating out the rotation spoils the 
linearity of the Hamiltonian in the lapse. To leading or¬ 
der in a, the rotation is ‘only’ linear in the lapse but to 
second order in a the non-linear dependence in the lapse 
becomes fatal. Integrating out the rotation, we find the 
Hamiltonian now carries a piece non-linear in the lapse 4 , 


H. 


total 


= Ho + HiN 


( 22 ) 


2 4 

OCpt 


64M|jTO 


10 [E]g 3 
O (R 3 /Mp^m 10 , a 3 ) . 


(e ij RiEj k Rk) 2 N 2 


As a result the Hamiltonian is actually genuinely non¬ 
linear in the lapse as soon as the non-minimal coupling 
is considered and the BD constraint is lost in the full 
theory. 

Extension to bi-gravity and multi-gravity. 

These arguments may easily be extended to multi¬ 
gravity theories. For instance in the case of bi-gravity, 
we must add to the set of auxiliary variables, the lapse 
and shift of the / metric, M,M l . The full set of auxil¬ 
iary variables are then the \D(D + 3) variables {y 1 } = 
{N^, M M , A b a }. Since bi-gravity caries one copy of dif- 
feomorphisms, we are guaranteed the existence of D first 
class constraints. This means that to ensure the absence 
of a BD ghost, we require that the Hamiltonian remains 
linear in the two lapses N° and M° and the shifts M l af¬ 
ter integrating out all the Lorentz Stiickelberg A b a . That 
this is not the case in the presence of matter coupled to 
the composite vielbein follows either by direct calculation 
or by recognizing that since massive gravity arises in the 
decoupling limit Mf —> 0, and in this limit M M decou¬ 
ples. Since this is a smooth limit, if the Hamiltonian was 
linear in the lapse N in bi-gravity then it would be so 
in the massive gravity limit which is a contradiction. In 
other words a BD ghost cannot reappear from a consis¬ 
tent decoupling limit. The same argument rules out all 
the multi-gravity extensions. 

Decoupling Limit. Although as we have emphasized, 
the unconstrained vielbein formulation is fundamentally 
different than the metric formulation, both have the same 
A decoupling limit. The vielbein derivation of the decou¬ 
pling limit including the vector degrees of freedom was 


4 Even though we have only solved the rotation to quadratic order 
in Ri it is sufficient to derive the Hamiltonian to quartic order 
in Ri because the Hamiltonian does not carry any piece in the 
rotation which is not also multiplied by it?. In other words, if we 
set Ri —> eRi and the rotation r — > e 2 r, then the Hamiltonian 
goes symbolically as H = 1 + e 2 ft? + e 4 (it| + rft?), so we only 
need to solve the rotation to quadratic order in Ri to know the 
Hamiltonian to quartic order in Ri. 
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given for massive gravity in [40, 41], bi-gravity in [42] 
and more recently in multi-gravity in [30]. To take this 
limit, we move out of unitary gauge and work with the 
reference vielbein in the form 


f a = (e A ) Q t> d<t> b , (23) 


where X a b are the Lorentz Stiickelberg fields and </>° the 
diffeomorphism Stiickelberg fields. We may for conve¬ 
nience choose e a to respect the symmetric vielbein con¬ 
dition. We are allowed to do this because on introducing 
the Lorentz Stiickelberg fields, we recover a local Lorentz 
symmetry which may then be used to set this condition 
on the vielbeins. With this condition, e a then encode the 
usual variables of the metric. Then the scalings which 
determine the decoupling limit, which as usual are deter¬ 
mined by canonically normalizing the quadratic fluctua¬ 
tions, are in general dimensions, 


e a = l a 


1 


( D - 2)/2 


M, 


pi 


1 


Xab — 


mM, 

1 


(D- 2)/2 


A a 


PI 


mM , 


(D-2)/2 

PI 


Xab 7 


(24) 

y\2+(.D-2)/2^ 7r >( 2 5) 
(26) 


were A( 2+ ( d_2 )/ 2 ) = m 2 Mpf 2 ^ 2 . In particular, in this 
limit Mpi —> oo for fixed a, the composite vielbein is 
given by 


e“ + ocf; = ((1 + a)(5“ - an%) + O (-) (?7) 

\mMp, / 

where II^j, = A _ ( 2+ (' D_2 )/ 2 )<9 /i d„7r. Following [22] we as¬ 
sume that the matter which couples to this effective met¬ 
ric, represented by y, does not scale. In which case, in 
the limit Mp\ —> oo, for fixed A, the matter Lagrangian 
remains finite with the fields y minimally coupled to the 
composite metric 

dfu = ((1 + ot)r]^ - alV) 2 , (28) 

and crucially A a and X a b drop out of the matter La¬ 
grangian. Hence the Lorentz Stiickelberg fields enter into 
the DL precisely in the same way as they do in the case of 
minimal coupling as given in [40, 41], which means that 
they continue to impose the DL version of the symmetric 
vielbein condition, i.e. the equation that determines A a b 
in terms of A a . The Lorentz Stiickelberg fields may then 
be integrated out to generate the standard vector inter¬ 
actions [40, 41]. The decoupling limit of the composite 
vielbein formulation is then identical the DL of the com¬ 
posite metric formulation [22], and the latter was shown 
to be ghost-free in [22]. This ensures that the ghost is 
absent below and around the scale A. Furthermore it 
is consistent to take Vainshtein screened solutions with 
n~ 1 without the ghost being a problem. 


Breaking the symmetric vielbein condition in 
multi-gravity. Another interesting situation where 
the symmetric vielbein condition is broken and the equiv¬ 
alence between the metric and vielbein formulation is 
spoiled arises in the context of multi-gravity [4, 43] and 
specifically in the presence of ‘cycles of interactions’. The 
metric formulation of such theories was conjectured to 
exhibit a ghost in [4] and this was proven directly in the 
decoupling limit in [43]. The vielbein formulation on the 
other hand breaks the symmetric vielbein condition. As 
such the equivalence with the metric is broken, leaving 
open the possibility that the vielbein formulation could 
be better behaved. To better understand the implica¬ 
tions of breaking the symmetric vielbein condition, we 
consider a specific example in D = 3 dimensions (see 
also [44]), 

C = MiCp,p[e\\ + M 2 Cp,p[e 2 ] (29) 

+M P \m 2 e A BC (cief A ef A ef 

+c 2 e A A ef A f c + c 3 f A A f B A ef) , 

where for simplicity we consider the “massive gravity” 
limit of tri gravity where one of the vielbein acts as a 
reference metric f A . We are dealing with a loop as long 
as none of the coefficients Cj vanish. If C 1 C 2 C 3 = 0 we 
are then dealing with a line of interaction rather than a 
cycle. 

Proceeding as before, we can work in unitary gauge 
for the reference vielbein f A = 6^ and decompose each 
of the two dynamical vielbein e\ and e 2 into a Lorentz 
transformation of an upper-triangular vielbein as in 
Eqns. (5-9). We can now proceed as earlier: one can 
use one set of shift equations to solve for one set of 
boosts. Then the remaining auxiliary variables, namely 
the remaining set of boosts and the two rotations can be 
integrated. 

In principle these integrations are quite involved but to 
go straight to the point we can perform a perturbation in 
the momentum constraints Ii ,. Then all the four boosts 
start at first order in Ri. To zeroth order in I?j, the 
Hamiltonian only includes the rotations and the lapses 
as auxiliary variables, 

H = iVi-R^o + N 2 R2,o + 2m 2 Mpi|ci[ei] + c 3 y/g 2 (30) 

+ Ni(ci + c 2 e iJ e a 6e li “e 2 j) + N 2 (c 3 [e 2 ] + c 2y fg{) j 
+0(Ri ), 

where we use the same notation as before, [e CT ] = e a \, 
and = e^ea^e^e^, and AI ct= i i2 are the two lapses 
(not to be confused with the shifts iV*). 

Since the boosts vanish to that order in Ri, the three 
symmetric vielbein conditions would read 

A = E 1 Z 1 E 1 2 ]S a b = 0 and n = r 2 = 0 . (31) 
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On the other hand, the equations of motion with respect 
to the rotation then leads to the modified symmetric viel- 
bein condition 

— C2Ni r )r 1 + r 2 '\/G'l 2 ^ (32) 

+ n + C2-/V1 (7^2 \/Gl2 + r 2 A)^j = 0 

- c 2 iVi 7 r2 ^7 ri „4 + riy/ Gi 2 ^ (33) 

+ r 2 ^C3-/V 2 [-B2] + C2N\{ r )r 1 \J G\2 — TiA)J = 0 , 

with \JG \2 = e u e a bE a °;E rj j. To that order, the sym¬ 
metric vielbein condition (31) is still a solution of the 
equation, but it is not the most general one. More gen¬ 
erally the two equations (32) and (33) can be solved for 
r\ and r 2 without specifying A. We can see immediately 
that if C 2 = 0, then the solution is simply r\ = r 2 = 0 
and we recover at least part of the symmetric vielbein 
condition. If ci = 0 or C 3 =0, one of the rotation still 
vanishes while the other one is uniquely determined in 
terms of the E a “ without invoking the lapses. Only in 
the case where C 1 C 2 C 3 ^ 0 and when A ^ 0 , the rotations 
depend non-trivially on the lapses. Integrating out the 
rotations we find that the Hamiltonian is no longer linear 
in the lapses and the theory excites the BD ghost 

H = NiRi t o + N 2 R 2.0 + 2 m 2 Mpi|c 3 V /g 2 + ci[i?i] 

+ N 1 (c 1 +c 2 VgT2) + N2(c2V9-i + c 3 [E2}) 

_ c 2 2 Nf(c 1 [E 1 }+c 3 N 2 [E 2 })A 2 _ 

2 {ciCzN 2 [Ei)[E 2 ) + C 2 V G 12 IV 1 (Cl[-El] + C 3 A ^2 [-^ 2 ])) 

+ O (A 3 ) } + O (. Ri ) . (34) 

The contributions arising from integrating out the boosts 
will all enter at least at first order in the momentum 
constraints Ri and cannot change the nature of the result 
obtained so far. We can see straight away that if we 
break the symmetric vielbein condition and 4/0, then 
as soon as C 1 C 2 C 3 ^ 0 , the second line leads to terms 
which are highly non linear in the lapses and thus spoil 
the constraints that should have been there to remove the 
BD ghost. The only possible way out is then to solve the 
equations of motion by imposing A = E x ,[ 1 ^, 2 ]^ = 
0. However this would overconstrain the system, and 
does not follow from any combination of equations of 
motion. As a result we confirm that cycles of interactions 
in multi-gravity break the symmetric vielbein condition 
and also lead to a BD ghost in the unconstrained vielbein 
language. This is consistent with the conclusions of [44]. 

Discussion. To summarize, the composite vielbein 
coupling is a ghost-free matter coupling in the decou¬ 
pling limit, where it is equivalent to the composite 
metric coupling considered in [ 22 ], however beyond the 
decoupling limit it propagates a BD ghost, independently 


on whether we consider it in the metric language or in 
the unconstrained vielbein one. Although we have not 
considered more general composite vielbein couplings 
here, any other composite metric or vielbein will lead to 
a ghost already in the decoupling limit 5 . 

These results may viewed positively in two ways, on 
the one hand there do exist ghost-free matter couplings 
in the decoupling limit, and that is a sufficient for 
validity of the EFT for energies E < A. On the other 
hand the absence of any non-minimal matter coupling, 
like the absence of non-minimal kinetic terms [45], is 
a remarkable testament to the power of field theory 
consistency in determining unique interacting theories. 

This analysis has also allowed us to study more 
general theories where the symmetric vielbein condition 
is broken and the vielbein formulation differs from 
the metric one. We argue that in that case the BD is 
also present in the vielbein formulation of the theory. 
Taken together with the results of [37] we are led to 
conclude that any consistent interacting theory for 
massive spin-two fields must have a natural expression 
in the metric language and that the vielbein formalism, 
although a very powerful tool from the calculational 
point of view and a very convenient formalism, does not 
give rise to new interactions. 
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